Phase Transition in Space: How Far Does a Symmetry Bend Before It Breaks? 
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We extend the theory of symmetry breaking dynamics in non-equihbrium second order phase 
transitions known as the Kibble-Zurek mechanism (KZM) to transitions where the change of phase 
occurs not in time, but in space. This can be due to a time-independent spatial variation of a field 
that imposes a phase with one symmetry to the left of where it attains critical value, while allowing 
spontaneous symmetry breaking to the right of that critical borderline. Topological defects need 
not form in such a situation. We show, however, that the size, in space, of the "scar" over which 
the order parameter adjusts as it "bends" interpolating between the phases with different symmetry 
follows from a KZM - like approach. As we illustrate on the example of a transverse quantum Ising 
model, in quantum phase transitions this spatial scale - the size of the scar - is directly reflected in 
the energy spectrum of the system: In particular, it determines the size of the energy gap. 
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INTRODUCTION 

Near the critical point of a second order phase tran- 
sition both the relaxation time (which determines "re- 
flexes" of the system) and the healing length (which sets 
the scale on which its order parameter "heals" in space, 
i.e., returns to its equilibrium value) diverge: 

r^Wkr, (1) 

c = vkr- (2) 

Above, e is a dimensionless parameter which measures 
the distance from the critical point. For instance, when 
the transition is caused by the change of the temperature 
T or by varying the parameter g of a Hamiltonian, e is, 
respectively, given by; 
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(3) 



where Tc and gc are the critical values. 

Divergences in t and ^ that appear near the critical 
point as a consequence of these two equations are often 
referred to as critical slowing down (for obvious reasons) 
and as critical opalescence (as the fluctuation on a scale 
^ that becomes large in the vicinity of a critical point 
cause, in some systems, variations of the optical proper- 
ties on scales ~ £,, which leads to opalescence). These 
divergences will play a crucial role in our discussion. 

Phase transitions are usually investigated as equilib- 
rium phenomena in homogeneous systems. But many in- 
triguing questions arise when a system is driven at a flnite 
pace from one phase to another. That this will inevitably 
happen in the cosmological setting was first pointed out 
by Kibble |l|, who noted that - because of relativistic 
causality - cosmological phase transitions in a variety 
of field theoretic models necessarily lead to formation of 
topological defects (such as monopoles or cosmic strings) 
and may have dramatic astrophysical consequences. 



It was later pointed out by one of us [2| that analogs 
of cosmological phase transitions can be studied in the 
laboratory, and that the equilibrium critical scalings - 
Eqs. ([T|) and ([2|)- can be used to predict various aspects 
of non-equilibrium dynamics of symmetry breaking, in- 
cluding the density of topological defects left behind by 
a non-equilibrium second order phase transition [2, y, Ij] . 

The resulting theory (known as "Kibble-Zurek mecha- 
nism" or KZM) uses critical scalings of the relaxation 
time and of the healing length to estimate the size ^ 
of the domains that choose the same broken symmetry 
[3, l3|. Owing to the universality of phase transitions, 
KZM can be applied on many energy scales, from the 
(cosmologically relevant) high energy settings (including 
experiments such as RHIC or LHC), all the way to Bose- 
Einstein condensates at ultra-low temperatures. 

In all of these situations a frequent prediction is that 
the process of symmetry breaking will lead to the forma- 
tion of topological defects. This is essentially inevitable 
when the broken symmetry phase (characterised by the 
homotopy group) permits their existence [ij, l2, y, Ij]. 
Following the transition, topological defects should ap- 
pear with the density of about one defect unit (e.g., 
one monopole or one ^-sized section of a string) per ^- 
sized domain. The value of ^ that results from this non- 
equilibrium process can be deduced from the equilibrium 
near-critical behaviour. Estimating this size is therefore 
essential. 

Here we first summarise the key ideas and equations 
that lead to such estimates. We shall then investigate the 
behaviour of the order parameter in an inhomogeneous 
system, where the change between phases is imposed by 
a slowly varying externally controlled parameter, such as 
the field in the quantum Ising model. We shall calcu- 
late how far does the transition region between the two 
phases persist - how much does a symmetry that prevails 
on one side of the critical line "bend" before it breaks, 
spontaneously, on the other side. 



DYNAMICS OF SYMMETRY BREAKING 

We consider a second order phase transition that is 
traversed at a finite rate set by the quench timescale tq , 



t 



(4) 



The system will be able to adjust its state adiabatically 
as long as the rate of change imposed from the outside is 
slow compared to its reaction time given by r, Eq. ([T]). 
This change from nearly adiabatic to approximately im- 
pulse behaviour will happen at an instant t when 



-(i") = ^, or rol-l— ^t 



(5) 



Thus, the state of the system's order parameter will in 
effect "freeze" at 



i^{roT!^n—=T. 



(6) 



Its evolution will restart only i after the critical point is 
passed. The instant t plays a key role in the establish- 
ment of the fluctuations which seed structures (such as 
topological defects) in the broken symmetry phase [5[- 
The key instant t corresponds to 






\TJ^ 



(7) 



which in turn sets the characteristic spatial scale given 
by the corresponding healing length 



To 



(8) 



This is the estimate of the size of regions that break sym- 
metry in a more or less coordinated manner [2, y] . Our 
derivation subverts equilibrium properties of the system 
- its scaling in the vicinity of the critical point ~ to pre- 
dict non-equilibrium consequences of the quench. The 
density of topological defects left behind in the wake of 
a phase transition is the best known (but not the only) 
example of such predictions of KZM. 

These predictions were tested, extended and refined 
with the help of numerical simulations [g, |7| , and verified 
in a variety of increasingly sophisticated and reliable ex- 
periments in liquid crystals [8|, |9| , supcrfluids [10|, |ll|, 13] , 
superconductors [Ij, [ij, llal , as well as other systems la |. 
The majority of the experiments to date are consistent 
with KZM. The case of superfluid ^He may be an ex- 
ception. There the initial reports of detection of KZM 
vortices [lO| were retracted [ll[ after it turned out that 
the observed copious vorticity was inadvertently induced 
by stirring. It is still not clear if *IIe behaviour is re- 
ally at odds with the KZM predictions re-evaluated in 
view of the refined numerical estimates [7[ • One problem 
is that the decay of vortex tangle generated by KZM 



may be faster than that generated by stirring: KZM 
leads to anticorrelations between neighbouring vortices, 
so phase ordering kinetics is expected to proceed faster 
than when vortices are correlated, as when they originate 
from turbulent fiows. An accessible recent summary of 
the experimental situation due to Kibble [l7| is a highly 
recommended and up-to-date complement to this brief 
discussion. Further experiments that may allow for a 
better control of various parameters would be obviously 
welcome. Gaseous Bose-Einstein condensates are one re- 
cent and very promising proving ground ISl. 191 that has 
already yielded some exciting results 20|, |2lJ, [2^ . 

Recently, KZM theory was applied to quantum phase 



transitions 
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There, the system is at T = 0, and the nature of its 
ground state changes discontinuously as a result of a con- 
tinuous change of some parameter of its Hamiltonian. At 
the critical value of that parameter the gap - energetic 
price of the lowest excitations, given by the difference 
between the energy of the ground and the first excited 
state - is at its minimum, and it disappears in the limit 
of an infinite system. In a finite system the quantum 
phase transition is then described by an avoided level 
crossing. The basic observation due to Damski [27| is 
that KZM can be used to analyse the avoided level cross- 
ing process, and that it yields results in excellent agree- 
ment with Landau-Zener theory [35|. As we shall see 
below, the analogy is based on the observation that, far 
away from the avoided level crossing, the quantum sys- 
tem usually starts its enforced evolution adiabatically, as 
was the case for second order phase transitions. How- 
ever, it can only react on the relaxation timescale given 
by the inverse of the gap. Therefore, near the avoided 
crossing - where the gap closes - its reficxes deteriorate, 
and the impulse approximation in the immediate vicinity 
of a quantum critical point is appropriate. Below we give 
an illustration of KZM in action in the exactly solvable 
and paradigmatic quantum Ising model. 

QUENCH IN A QUANTUM ISING MODEL 

According to Sachdev [36| the quantum Ising model is 
one of two prototypical examples of quantum phase tran- 
sitions. It represents a chain of spins with the Hamilto- 
nian 



H = -Jit) 
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(9) 



where af'^ are Pauli operators, W is the Ising coupling, 
and J{t) is due to the external (e.g., magnetic) field that 
attempts to align all spins with the a;-axis. The phase 
transition from the paramagnetic state where all the 
spins are aligned with x by the strong external field (e.g., 
I ^■, —>,... , — >)) to the low-field ferromagnetic and de- 
generate (in the large TV limit) ground state that "lives" 



in the Hilbert subspace spanned by the broken symme- 
try basis I t, t, • ■ ■ , T) and | i, i, ■ • ■ i) takes place when 
J{t) = W. Therefore, as suggested by Eqs. ([2]and|l]), a 
key role is played by the relative coupling given by 



S{t) = J{t)/W - 1 = t/TQ 



(10) 



and we assume that it depends linearly on t. 

All the relevant properties depend on the size of the 
gap A between the ground state and the first excited 
state. In an infinite system the gap is given by: 



A = 2|J(i)-W^| ^2W\e{t)\ 



(11) 



The spectrum (and the gap) of a finite quantum Ising 
system are illustrated in Fig. 1. The gap sets the en- 
ergy scale and is reflected in the relaxation time and the 
healing length 

T = h/A = h/2W\e{t)\=To/\e{t)\, (12) 

e = 2Wa/A{t) = a/\e{t)\ = ^o/\e{t)\ , (13) 

the latter given by the product of r and the speed of 
sound c = 2Wa/h (see [36|), where a is the distance 
between spins. The divergence of r and ^ is the critical 
slowing down and the analogue of critical opalescence, 
respectively. 

The instants ±t when the behaviour switches from adi- 
abatic to impulse happen when the reaction time of the 
system, Eq. ^T2\i . is the same as the timescale on which 
its Hamiltonian is altered (given by e{t)/e{t) = i), i.e. 



Tii)^ To/\e{i)\ ^ e{i)/iii)^i, (14) 



and thus 



i=^T^, = ^TQh/2W. (15) 

As before, for t < — i, the state of the system will con- 
tinue to adjust adiabatically to changes imposed by the 
decreasing J{t). However, at t = —t before the critical 
point the evolution will cease, and it will re-start only at 
t = +t after the transition, with the initial state similar 
to the one "frozen out" at — i. 

Using the relative coupling e associated with t we get 



e = e(i) = t/TQ = Jto/tq 



(16) 



( = Ca/e = Ca^jTQ/To = a^2WTQ/h. (17) 

Following KZM, wc now predict appearance of 0(1) de- 
fects per ^. Their density should be approximately 



fiKZAi ^ a/C = \ f^f^WTQ 



(18) 



per spin. This is only an estimate: Simulations of classi- 
cal second order phase transitions yield defect densities 
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FIG. 1: Energies of lowest excitations of the Ising chain for 
A*' = 20 with open (rather than periodic) boundary condi- 
tions. The ground and the first accessible excited state that 
define the gap are marked with a thicker line. 



that scale with tq in accordance with this reasoning, but 
are lower than the relevant power of ^ of Eq. ^T7\ so that 
a "unit of defect" is often separated by ^--^10-15 £,; see e.g. 

There is a fundamental difference between quantum 
and thermodynamic phase transitions. In the thermo- 
dynamic case 'real' fluctuations exist above the critical 
point. They can initiate the symmetry breaking pro- 
cess " in effect, choose how symmetry breaks in domains 
that appear after the transition. In the case of quantum 
phase transition there are only 'quantum fluctuations', 
but they are virtual, so one cannot be certain that they 
will have an analogous effect on the post-transition state. 
Therefore, an explicitly quantum approach to the quan- 
tum Ising model is needed. 

As we have already noted, a key feature of a quantum 
phase transition is the gap A. Actually, as can be seen in 
Fig. [U the relevant gap (i.e. the gap between the ground 
and the first accessible state) is not the symmetric A but 
rather 2A on the approaching side. Nevertheless, in the 
quantum Ising model the gap disappears at the critical 
point when the system is infinite. When N < oo, this 
critical gap becomes small, but docs not disappear, (see 
Fig. [T]). This is important, and allows us to propose a 
purely quantum approach: Instead of density of defects 
in an infinite system wc can compute (as a function of 
quench timescale tq) the length of the spin chain that, 
with probability p of about a half, can remain defect - 
free (in a ground state) after the quench. An excited 
state would (most likely) contain just a single excitation 
(e.g., a kink). Therefore, the inverse of the length that 
limits the excitation probability to approximately '^ 0.5 
should correspond to about half of the kink density. 

The two lowest accessible levels of H in the vicinity of 
the critical point (Fig.[T]) exhibit an avoided level crossing 



and we can calculate the excitation probability of the 
system driven through this avoided crossing using the 
Landau-Zener formula (LZF) [3a |. 



P' 



(19) 



Here A is the minimum gap between the two levels and 
V is the velocity with which the transition is imposed on 
the system given by u = A far away from the avoided 
crossing. The KZM approach was shown by Damski to 
provide an excellent approximation to LZF [27| (see also 
|28| for extensions). Using LZF wc compute the size N 
of the spin chain that will probably remain in the ground 
state in course of the quench with probability p ~ 0.5. 
Equation p9|) translates into a condition for the rate of 
quench that produces a kink with probability p, 



\v\< 



7rA2 



2;i|lnp| 



(20) 



Using w = |A| = 2J(i) = ^W/tq [see Eq. ^] and 
A = AttW/N for the gap upon "closest approach" we get 



lAI 



2W Tr{ATrW/N)^ 



2h\\np\ 



(21) 



This relates the size TV of a chain that will remain defect- 
free with the probability 1-p to the quench rate: 



A^< 2n, 



IttWtq 
h\ \np\ 



(22) 



This LZF estimate is (surprisingly) accurate for p < 0.5 
even though there are many levels in the spectrum of 
the quantum Ising model. Such accuracy was surpris- 
ing when this analysis was first presented |23| . Dziar- 
maga [25| has soon after demonstrated that this is no 
accident - in effect, the phase transition in the quantum 
Ising model can be decomposed into a collection of inde- 
pendent avoided level crossings. 

Wc can now directly compare KZM, Eq. p8)) . and LZF 
predictions for defect density: 
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FIG. 2: (a) Number of kinks per spin in the quantum Ising 
model after a quench starting in the ground state at J = 5W 
and ending at J = 0, plotted as a function of the quench rate 
tq/tq = hv/AW'^ for iV = 50,60,70,80, 90,100 (solid lines; 
bottom to top). The scaling ukzm ~ \/to/tq predicted by 
KZM is consistent with the simulations (see [33| for details 
of the numerical method) . Agreement improves with the size 
of the system: for 100 spins a fit gives n ~ "'"o"'^* (dashed 
line). As in the classical case [g] Eq. (|18|l is an overesti- 
mate; the best fit is n ~ O.lQfiKZM ■ (b) Total number of 
kinks for A^ = 50, 60, 70, 80, 90, 100 spins (dashed lines, bot- 
tom to top) after a quench as a function of the quench rate 
T()/t"q = hv/4W^. Both the scaling ukzm ~ l/y/rQ pre- 
dicted by KZM (solid, straight lines), Eq. HHJ, and the LZF 
estimate p (solid, bend lines) when less than one kink is ex- 
pected are valid. The straight lines are linear fits in the range 
(0.025,0.25) yielding slopes between 0.66 and 0.58. Numerical 
data include these used in (a) but we now go beyond the ex- 
pected range of vaUdity of KZM. For sufficiently slow quenches 
LZF provides reliable predictions. Very fast quenches are "aU 
impulse" , levelling off of the expected number of kinks, as is 
indeed seen. 
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(23) 



The two estimates exhibit the same scaling with the 
quench rate and with the parameters of iJ, Eq. ^. 
LZF predicts fewer defects than "raw KZM" {fiLZF — 
0.14 X fiKZM when p is set - somewhat arbitrarily - to 
0.5). This is not a big surprise - as seen in the numerical 
simulations, confirmed by the experiments and verified 
analytically in specific models, Eqs. (fTT]) and (fTS]) pro- 
vide correct scalings, but tend to overestimate densities 
(see e.g. [y, lZ[). Fig. [5] indicates that this conclusion 
holds also for the quantum Ising model. 



Figure [5] also shows that the kink density scales ap- 
proximately as ^ \/ y/TQ, Eq. p^ . in the region of the 
validity of KZM. i.e. for e less than 1 (so that quench 
is quasi - adiabatic at the beginning and at the end, but 
impulse near the critical point, i.e., when at least one de- 
fect is expected). The prefactor which is approximately 
0.16 [0.12 if the steeper slope on the approach in Fig. [T] 
is taken] is not far from the previous experience (6|. For 
very slow quenches {i > h/A, or tq > (^)^to), i.e. if 
the system is nowhere convincingly 'impulse', LZF is sur- 
prisingly accurate. We conclude that the two approaches 
work well in complementary regimes of quench rates, and 



predict the same scaling of the size of broken symmetry 
domains with quench time. 

The importance of the behaviour of the gap (and, in 
particular, of the smallest gap A) for the quantum es- 
timate of defect density is - in addition to the broader 
theme of connections between KZM and LZF - our prin- 
cipal general conclusion at this point. We shall return 
to discuss the gap in the quantum Ising model and util- 
ity of extensions of KZM and LZF to more complicated 
behaviours of the near-critical gap below. 

STATICS OF SYMMETRY BREAKING 

Consider now a situation where the all-important di- 
mensionless e is time-independent, but depends on the 
location in space instead. We suppose - in analogy with 
Eqs. (HI [TO]) ~ that in the vicinity of the critical point e{x) 
is approximately linear: 



from the critical point. The distance x corresponds to 
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(24) 



Clearly, sufficiently far from x ~ the order parameter 
will settle into an equilibrium state corresponding to the 
local value of e. We are however left with an interesting 
question: How far - and how - does the influence of the 
critical point propagate? 

There is an intriguing analogy between the freeze-out 
in time we have discussed in the preceding section and 
the freeze-out in space we are led to consider here: When 
e changes in space slowly compared to the local healing 
length given by ^ = ^o/kT; Eq. ^, the order param- 
eter will be able to adjust in space to the changes of e 
- i.e., there should be a well-defined local value of the 
order parameter, local healing length ^, and local relax- 
ation time r, etc. However, very close to a; = the 
system cannot "heal" fast enough: The healing length 
becomes large, much larger than x. Consequently, the 
critical opalescence - the divergence of ^ = ^o/leT - will 
open up a "scar" in the order parameter that does not 
properly heal. Our aim now is to describe the conse- 
quences. To find out the size of the scar we write down 
the spatial analogue of Eq. ^ , 



ai) 



e{x) 



for the "adiabatic-impulse" borderline point, x. Thus, 



(25) 



Coif I 



(26) 



Consequently, and in accord with previous discussion, we 
conclude that the order parameter will in effect "freeze" 
(or, to put it differently, that the scar opened up by the 
transition through the critical region will heal) at a dis- 
tance 



(f)^- 
^Q 



(28) 



We are really done - there is just a spatial x, rather than 
both X and t, as before. But let us (as a consistency 
check) repeat and take a few more steps, and use this 
estimate of e to calculate the characteristic spatial scale 
given by the corresponding healing length. 



C = Co( 
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(29) 



We know this already [sec Eq. (p7)) ] . This is now the es- 
timate of the size of the scar - the size of the region that 
is "still thinking" about how to break symmetry. Earlier, 
our derivation of the freezeout time and of the resulting 
analysis of the dynamics of symmetry breaking subverted 
equilibrium properties of the system - its scaling in the 
vicinity of the critical point - to predict non-equilibrium 
consequences of the quench. We have now repeated this 
subversive strategy in a new setting, by using the healing 
length (that is, strictly speaking, defined in a homoge- 
neous system) to find out consequences of inhomogeneity 
that must be there in the vicinity of the critical point if 
the phase transition takes place in space. The structure 
of the scar in the order parameter that connects the two 
phases is now of interest per se. 

There are two specific values of ly that are often en- 
countered: As we have seen in the preceding section, for 
quantum Ising model ly ~ 1. Therefore, 



x^ VCoAq = C 



(30) 



and the size of the scar is simply a geometric average 
of the two relevant lengths in the problem. Wc shall 
see below that, as a consequence, in the quantum Ising 
model the presence of such a scar widens a critical gap. 
In a mean field theory we have i^ = 1/2. Consequently, 



(eo'AQ)^ = e 



(31) 



(CoA^) — 
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(27) 



The thickness of the boundary between the two layers 
is a compromise that involves the characteristic healing 
length of the order parameter on one hand, and the typ- 
ical scale of the externally imposed inhomogeneity of the 
Hamiltonian on the other. The extent to which each of 
them has its say depends on the scaling exponent v. The 
most obvious test of these predictions would then explore 
variations of the order parameter in the vicinity of an ex- 
ternally imposed inhomogeneity. 

OPENING A GAP WITH A WEDGE OF A FIELD 

To illustrate our considerations with a concrete exam- 
ple we return to the quantum Ising model we have already 
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FIG. 3: Illustration of the spatially dependent field J(x) in 
the transverse quantum Ising model we consider. The param- 
eter J(x) corresponds to J(f) from Hamiltonian ((Qjl but now 
depends on position, J{t) -^ ,J(x). The spins are located at 
positions xi = (I — l)a with I = 1, . . . ,N. 



employed in the previous section. The Hamiltonian we 
shall use now has the same form as before, Eq. ^, but 
the field J, which before was the same for every spin, now 
has a form illustrated in Fig. [31 As a result, sufficiently 
far to the left of x the system will be in a paramagnetic 
phase, while to the right a ferromagnetic phase should be 
dominant. 

Our task is to characterise the state of such a system. 
In a homogeneous quantum Ising spin system the corre- 
lation function given by 



C(^-) 



l-k'^l 



(32) 



is a well defined object, whose behaviour is thoroughly 
explored [36[. Indeed, the decay of such correlation func- 
tions is used to define the coherence length (which is in 
effect often equal to the healing length we have employed 
earlier). In a homogeneous system ((Tf_j,crf) obviously 
does not depend on the "reference spin" I: Each spin 
"lives" in an identical neighbourhood, so C{k) is transla- 
tionally invariant with respect to I. But when a system 
is inhomogeneous - for example when a phase transition 
occurs at some point within the system - the correlation 
function will obviously depend on where the reference 
spin I is. We therefore consider a conditional correlation 
function, 



mi) 



'i-k"i ^ 



(33) 



It depends explicitly on the location of the reference spin 
I. Formally, C,{k) is given by ^J.{k\l) averaged over all /. 
In a homogeneous case ^i{k\l) is independent of ?, so this 
averaging is trivial. On the other hand, in an inhomo- 
geneous system we are investigating here, such averaging 
would be counterproductive: The averaging would ob- 
scure precisely the signature we are looking for, i.e., the 
imprint made by the inhomogeneity of the Hamiltonian 
on the state of the system. 

In Fig. [4^ we show n{k\l) for the fixed reference spin 
/ located in the ferromagnetic region. The behaviour of 
pL{k\l) seen there is a good illustration of our previous 
discussion: For spins located in the ferromagnetic region 
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FIG. 4: (a) Conditional correlation function, /i(fc|Z), for A'^ = 
200, x/a = 120, aa = 0.05 and I = 160. The reference 
spin located ai x = xi is in the ferromagnetic region, (b) 
Width of the transition region of the conditional correlation 
function i-i{k\l) versus a for A'^ — 200, x/a — 120 and I = 
160. The dashed line is a fit to the data leading to La = 
2.34(Qa)-'^-54 = 2.3A{\Q/a)°-''\ 



fj.{k\l) w 1. However, when k assumes values such that 
I — k falls within the paramagnetic region, the correlation 
function vanishes, fJ-{k\l) « 0. The transition between 
these two extremes is our focus of interest. We can char- 
acterise the size of the "scar" by quantifying the width 
of the transition regime. The most obvious way to do 
this is through the inverse of the slope of fJ.{k\l) where 
/i(fc|/) = i The resulting width of the transition region 
is plotted in Fig. HJd. Evidently, the behaviour predicted 
by the KZM-like discussion for the quantum Ising model, 
Eq. ([501) . is refiected in the correlations between spins. 
Such correlation functions should be experimentally ac- 
cessible, and may be a good way to test our predictions. 
We note that some of the effects of the inability of the 
order parameter to adjust to the inhomogeneous varia- 
tions of the Hamiltonian may become apparent even if the 
system does not cross "all the way" into the other phase. 
For instance, just coming close to the critical point may 
lead to the structures with sizes that can be estimated us- 
ing the above approach. Similarly, when the same phase 
is separated by a narrow strip when the parameter that 
controls its phase crosses the critical point (i.e., a strip 
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FIG. 5: Energies of lowest excitations for A'^ — 50, a — 0.05/a 
versus x. The ground and the first accessible excited state 
that define the gap are marked with a thicker line. 



in space within which e has a different sign than outside 
it), an extension of the above discussion should be ap- 
pUcable. There should be therefore lots of opportunities 
to experimentally test the above equations, and plenty 
of variations of the homogeneity - inhomogeneity theme. 
One such obvious variation involves adopting a nonlinear 
spatial dependence of e on x - i.e., a dependence that is 
different from Eq. (j24p . But there are clearly many more. 
Wc shall not attempt to enumerate them here. 

Let us instead point out a less obvious consequence of 
a spatial inhomogeneity that will be especially important 
in quantum phase transitions: Its effect on the gap. The 
spectrum of an inhomogeneous Ising system as a func- 
tion of the location of the point x - the place where J{x) 
becomes 0, see Fig. |3]- is shown in Fig. [5l When this 
figure is compared with Fig. [U the spectrum of a homo- 
geneous system, two striking differences in the behaviour 
of the gap emerge, (i) Instead of a sharp "corner" where 
J/W = 1, there is now a plateau which extends over a 
range of values of x. (ii) Moreover, this plateau is lifted 
above where it was (at ^^ AirW/N) for a chain of A^ spins. 
In short, the sharp gap minimum we had before becomes 
now a wider, extended "bottleneck" . The spatial struc- 
ture induced by the inhomogeneity in the state of the 
system should have an effect on their energetics. Indeed, 
one could venture a guess (based on "general principles" ) 
that the size of the gap - which is the property of the 
whole inhomogeneous spin chain - will change with the 
size of the scar, x. One could moreover speculate (and 
here the "general principles" are becoming too general 
for comfort) that the size of the gap should be inversely 
proportional to (some power of) x. Wc shall now pose 
this very same question more precisely, and arrive at the 
estimate of the basic parameters of the "bottleneck" more 
directly. 

The inhomogeneous Ising system we are dealing with 




FIG. 6: Illustration of inequality (|34|) which defines a critical 
region of length Xa — zi ~ zq = 2A/aW. 



can be typically divided into three regions: A ferromag- 
netic domain (where the Ising coupling W dominates), 
a paramagnetic domain (where J{x) is large) and what's 
in between these two - the near-critical domain, where 
\J{x) — W\ is small. The energetic price of excitations 
in either ferromagnetic or paramagnetic domains is large 
compared to the size of the gap, which (we remind the 
reader) in a near-critical chain of L spins is given by 
Al ~ AttW/L. 

Our aim is to compute the size of the bottleneck gap 
in the inhomogeneous case of Fig. [31 We shall do that 
by a self-consistent estimate suggested by the above di- 
vision of the whole inhomogeneous system into three do- 
mains. We first note that the energetic price of the ex- 
citations outside the near-critical regime is prohibitively 
large whenever these regimes are well defined, so such 
excitations will not be relevant for the calculation of the 
size of the gap (which is defined as the "price" of the least 
energetically expensive excitation). On the other hand, 
the near-critical domain should be able to support col- 
lective excitations that - when the size of that domain 
is L - should be able to support collective eigenstates 
with energies given by A^ ~ AttW/L. Thus, all we have 
to do is to find that size L. One might be at this stage 
tempted to venture a guess that L ~ a:, but let us proceed 
cautiously, and follow a line of reasoning based on ener- 
getics we have outlined above. To this end wc propose a 
self-consistency condition. 



\J{x)~W\ < A, 



(34) 



which in effect defines the near-critical domain. A plot 
that illustrates this inequality is shown in Fig. [B) 

At this point the relevant A is not known. But (as 
\J{x) — W\ is small in the near-critical domain) it is nat- 
ural to employ Al ~ 4ttW/L, where now L stands for 
the number of spins that contribute to the collective state 
that defines the gap, i.e., the spins that fit inside the near- 
critical domain. Using J{x)/W = —a{x — x) of Fig. [31 we 
solve the above equation to obtain the extent, in space, 
of the near-critical domain that defines the size of the 
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FIG. 7: The size of the smallest gap between the ground 
state and the first accessible excited state - the level of the 
plateau in Fig. [5]- as a function of the slope a = 1/Aq for 
A'' = 50 spins. The dotted line is a fit to the data between 
aa = 0.01 and aa = 0.2 leading to l^^i/W = 2.72(aa)°''^ 
The horizontal, dashed line corresponds to the minimum gap 
of a homogeneous spin chain. 



gap, 



Xa = v/87ra/"~ 5.01a 



(35) 



or, in number of spins (rather than in distance). 



La = X^/a = yJ^-K/aa w 5.01/ 



(36) 



This immediately yields an estimate of the gap size, 
A^ = A-kW/La ~ 2My/mW. 



(37) 



With the help of these results we now understand the 
basic structure of the eigenenergies, and, in particular, 
of the gap seen in Fig. \5[ Figure [7] shows the size of 
the gap as a function of the slant a = 1/Aq. The square 
root dependence predicted by Eq. ([57)) is evident. Indeed, 
even the prefactor obtained through our simple estimate 
above is close to the one that is obtained from the numer- 
ical experiment. There are obvious limits to the validity 
of our derivation that are also in evidence with Fig. [T] 
When a is so small that the system as a whole is approx- 
imately homogeneous, the near-critical domain extends 

over the whole system, \/ ^ ~ -^A is comparable to, 
or larger than N, and the asymptotic size of the gap 
^ AttW/N for the homogeneous system is attained. On 
the other hand, when a is so large that the size of the 
scar becomes comparable with a, the spacing between 
the spins, the near-critical region disappears altogether, 
and the size of the gap is set by 2W. 

DISCUSSION 

We have presented an approach to symmetry breaking 
phase transitions that occur in space. As in the original 



KZM (which employs equilibrium near-critical scalings 
of the relaxation time r and of the healing length ^), it 
rests on the near-critical behaviour of £,. While the heal- 
ing length is defined in a homogeneous system, we have 
employed its variation with e (dimensionless parameter 
that measures the distance from the critical point) to 
predict what happens near a critical point that occurs 
in space in an obviously inhomogeneous system. The re- 
sulting theory yields the size of the scar of the transition 
region between the two phases. 

We have applied these results to the quantum Ising 
model. There the gradient of the external field leads 
to a gradual re-alignment of the spins. The border be- 
tween the paramagnetic and ferromagnetic phases has a 
structure that can be characterised by the correlations 
between spins, which change between the two patterns 
of alignment. We have seen that this change is grad- 
ual, and happens (as was predicted by our theory) on 
the scale inversely proportional to the square root of the 
gradient of the field. The simple theory leading to Eq. 
([29| results in correct scaling behaviour. As expected, 
a more specific calculation based on the quantum Ising 
model leads to a different prefactor, Eqs. ((35l [36)) . In 
effect, we now have two different measures of the size of 
the scar, one based on the correlations. Fig. [H and the 
other on the size of the gap. We do not know which one 
of them (if any) is "correct" , but we point out that they 
can differ. 

Existence of the "transition scar" has a dramatic ef- 
fect on the eigenspectrum of the system. In particular, 
instead of the simple avoided level crossing where the 
minimum gap helps set the probability that the system 
will remain in its ground state we are now dealing with a 
an extended "bottleneck gap" . Parameters of this bottle- 
neck gap are (including a prefactor, somewhat in contrast 
with the case of the scar size estimates) in a surprisingly 
good accord with the theory. 

One is now tempted to consider quantum phase transi- 
tions that occur both in space and in time. For thermo- 
dynamic phase transitions this situation can be treated 
by the KZM-based approach of Kibble and Volovik [38| . 
Their conclusion was, in effect, that when the critical line 
moves through space with velocities in excess oi v = ^/f, 
the quench that leads to spontaneous symmetry break- 
ing will proceed as if it were homogeneous, and KZM 
estimates of defect densities relevant for homogeneous 
quenches will apply. On the other hand, in the limit of 
slowly propagating phase fronts defect production will be 
suppressed, as the symmetry breaking choices made by 
the order parameter will be simply propagated with the 
phase front. This conclusion seems to be only weakly 
influenced by the spatial gradient of e, and is now sup- 
ported by numerical studies (see e.g. |39j). 

Our preliminary numerical simulations of the same 
problem in the quantum case indicate that there the sit- 
uation is more complicated. Part of the problem can 
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FIG. 8: (a) Energy of the first accessible state for N = 80, 
a — 0.05/a (solid line). The slant is moving with a constant 
velocity a/rg across the spin chain, i.e. x{t) = xq — at/rQ. 
xo is chosen such that the flat part of the spectrum begins at 
i = 0. The dashed lines are linear fits to the solid line in order 
to get parameters for a Landau-Zener type two-level Ifamilto- 
nian the excited state of which is given by the dotted line, (b) 
Upper line: Solution of a Landau-Zener type problem with an 
energy spectrum given by the dotted line in (a). The func- 
tion oscillates rapidly around the value 7.2 x 10~'^ (H/Wtq)^ 
(dotted line). The lower line shows the excitation probability 
for A'' — 80, a = 0.05/a. The dashed line is a fit leading to 
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be traced to the fact that the dynamics of a quantum 
phase transition is strictly reversible. Thus, the energy 
deposited in the order parameter stays in the order pa- 
rameter, and - e.g., in the quantum Ising model - has 
"nowhere to go" except into the kinks. As a result, the 
rate of creation of kinks now depends on the spatial gradi- 
ent of e, and does not become suppressed as dramatically 
as before when the front velocity falls below v. Moreover, 
in quantum phase transitions, the limit when less than a 
single kink is expected in the whole system is explicitly 
quantum, and (as we have seen earlier) different than (al- 
though compatible with) the KZM scaling. However, in 
an inhomogeneous case the simple Landau-Zener formula 
(suitable for gaps of Fig. 1) is obviously inapphcablc. We 
have tried out its natural generalisation (which approx- 
imates the spectrum of Fig. [5] with two obvious avoided 
level crossings connected with a plateau (see Fig. |5^). 



The fit is imperfect. The resulting prediction of the prob- 
ability of a kink. Fig. [SJd, is also not very satisfying. But 
this calculation (while it predicts spurious oscillations, 
and overestimates the rate by about an order of magni- 
tude) does capture one essential feature: The dependence 
of the probability p on the quench rate 1/tq is now no 
longer exponential, as it was in the original LZF, but it 
becomes quadratic. Such quadratic dependence appears 
automatically, if rather surprisingly - given the expo- 
nential nature of LZF, Eq. P^ - in the avoided level 
crossing problem when the transition starts at the place 
of the nearest approach of the two levels 28|, |30| , rather 



than far away from the avoided crossing, as is usually 
assumed [35|. The behaviour we observe may be due 
to the asymmetry of the gap. In effect, as the slanted 
potential of Fig. [3] traverses the spin chain, the system 
will generally approach the plateau of the bottleneck gap 
with a slope that is quite different from the slope on the 
other end of the gap, where it exits. Therefore, recent 
study [2^ - showing that when avoided level crossing is 
traversed starting at the minimum of the gap, the proba- 
bility of excitation p on the quench rate 1/tq is no longer 
exponential but that it becomes quadratic - becomes rel- 
evant. For instance, in the case of the asymmetric gap 
of Fig. [HK, one can imagine that the system is delivered 
essentially in its ground state to the bottleneck of the 
gap. While x moves along the fiat part of the gap, the 
Hamiltonian of the system is unchanged, so no transi- 
tions happen. However, as the gap opens up, now with a 
steep slope, one is really starting an avoided level cross- 
ing transition from a system that is still in its ground 
state near the end of the plateau. Clearly, conclusions of 
Refs. ^^ apply 

While these preliminary insights are encouraging, more 
work is needed to gain a complete picture. In addition 
to the obvious intrinsic interest of this problem there are 
applications that may benefit from its thorough under- 
standing. Let us mention one, slightly speculative: Adia- 
batic quantum computing rests on the idea that a known 
ground state of a simple Hamiltonian can be adiabati- 
cally transformed into an initially unknown ground state 
that solves a problem encoded into the structure of an- 
other Hamiltonian. The rate at which this transition 
can be accomplished is limited by the size of the gap - 
faster transitions will result in errors appearing at the 
rate given by LZF, Eq. p^ . This will severely limit the 
speed of the adiabatic computation - its duration will 
increase exponentially with the inverse of the size of the 
gap. But we have seen above that a gap can be widened, 
and the rate of transitions can be brought down from 
exponential to quadratic in 1/tq. It remains to be seen 
whether (and how) these changes can be used to suggest 
improvements in adiabatic quantum computing, but the 
preliminary results we have reached here are certainly 
suggestive. 
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